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ABSTRACT 

The paucity of old isolated accreting neutron stars in ROSAT observations is used to derive 
a lower limit on the mean velocity of neutron stars at birth. The secular evolution of the 
population is simulated following the paths of a statistical sample of stars for different values 
of the initial kick velocity, drawn from an isotropic Gaussian distribution with mean velocity 
< (y) < 550 kms~^. The spin-down, induced by dipole losses and the interaction with the 
ambient medium, is tracked together with the dynamical evolution in the Galactic potential, 
allowing for the determination of the fraction of stars which are, at present, in each of the four 
possible stages: Ejector, Propeller, Accretor, and Georotator. Taking from the ROSAT All Sky 
Survey an upper limit of ~ 10 accreting neutron stars within ~ 140 pc from the Sun, we infer 
a lower bound for the mean kick velocity, (y) ^ 200 — 300 kms~^, corresponding to a velocity 
dispersion ay ^ 125 — 190 km s~^. The same conclusion is reached for both a constant magnetic 
field {B ~ 10^^ G) and a magnetic field decaying exponentially with a timescale ~ 10^ yr. Such 
high velocities are consistent with those derived from radio pulsar observations. Present results, 
moreover, constrain the fraction of low velocity stars, which could have escaped pulsar statistics, 
to less than 1%. 

Subject headings: Accretion, accretion disks — stars: kinematics — stars: neutron — stars: 
magnetic field — X-rays: stars 



1. INTRODUCTION 

Isolated neutron stars (NSs) are expected to be as 
many as 10*^-10^, ~ 1% of the total stellar content of 
the Galaxy. Young NSs, active as pulsars, comprise 
only a tiny fraction (~ 10~^ — 10~^) of the entire pop- 
ulation, and about 1,000 have been detected in radio 
surveys. It is thus an observational fact that most 
of the NSs remain undetected as yet. Despite inten- 
sive searches at all wavelengths, only a few (putative) 
isolated NSs which are not radio pulsars (or soft 7 
repeaters) have been recently discovered in the X- 
rays with ROSAT (jWaTter, Wolk fc Neiihauser 199^ ; 
Haberl et al. 1996|; Paberl, Motch fc PietschT998| ; 



part has been found in the X-ray error box of RX 



Schwope et al. 1999| ; [Motch, et al. 1999| ). All these 
sources emit a thermal spectrum at ~ 100 eV and 
the derived column densities place them at relatively 
close distances (at most a few hundred parsecs) im- 
plying luminosities not in excess of 10^^ ergs~^. For 
RXJ1856-3754, an extremely faint optical counter- 
part {my ~ 26) have been firmly established by now 
( [Walter fc Matthews 1997| ) while a possible counter- 



J072 0-3125 by [Kulkarni k van Kerkwijk 1998[ (see 
also Motch &: Haberl 1998 ). For the other remain- 
ing condidates, plate searches have only placed limits 
down to my ^ 24. 

Although the extreme X-ray to optical flux ra- 
tio (> 10'^) makes the NS option rather robust, the 
exact nature of their emission is still controversial. 
Up to now, two main possibilities have been sug- 
gested, either relatively young NSs radiating away 
their residual internal energy or much aged NSs ac- 
creting the interstellar gas. Unmagnetized, cooling 
models can reproduce the full spectral energy distri- 
bution (SED), including the optical excess observed 
in the best studied object RX J1856-3754, if the 
surface temperature depends on the latitude or for 



Fe/Si-ash atmospheres (Walter & Matthews 1997; 
[Walter fc An 1998 ; note, however, that spectra com- 
puted by Pavlov, et al. 199^ for similar compositions 
fail to reproduce the observed SED). Accreting, pure 
H, atmospheres around non-magnetic NSs produce a 
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hard X-ray tail ( Zampieri et al. 1995| ) that is not ob- 
served in these sources. However, as in cooling mod- 
els, the hardening decreases for increasing B and the 
essential features of the SED may have a quite natu- 
ral explanation in terms of magnetized, H, accretion 
models (|Zane, Turolla fc Treves 1999| ). 

Cooling NSs have short lifetimes (~ 10^ yr; Page| 
1998|}| and thus might be relatively rare objects. On 
the other hand, at least if the Bondi-Hoyle scenario 
applies, accretion can be severely reduced when the 
NS velocity (relative to the interstellar medium) ex- 
ceeds 40 kms~^ so the process of accretion itself 
might be unable to produce the typical luminosities 
inferred from ROSAT data. When available, proper 
motion measurements for some of the five isolated NS 
candidates detected so far will prove decisive in as- 
sessing the true nature of these sources. Meanwhile, 
we feel that a more thorough analysis of the statisti- 
cal properties of NSs is of interest and can be useful 
in providing indirect evidence in favor of or against 
the accretion scenario. 



ever, little attention was paid to the NSs magneto- 
rotational evolution. Only recently, this issue was 



discussed in some detail by Livio, Xu fc Frank (1998) 



and |Colpi et al. (1998)|, who found that, for a given 



velocity distribution, the number of accreting sources 
depends strongly on the way the magnetic field de- 
cays. 

Goal of this investigation is to consider these two 
issues simultaneously, coupling the dynamical and 
the magneto-rotational evolution for the isolated NS 
population. In particular, we explore the effects in- 
duced on the current census of NSs by varying the 
mean value of the kick velocity. In earlier investi- 
gations (pTreves fc Colpi 1991| ; |Blaes fc Madau 1993| ; 
Zane et al. 1995|), the velocity distribution at birth 



was assumed either Gaussian with a mean veloc- 
ity modulus (inferred on the basis of existent data, 
Narayan fc Ostriker 1990| ) of ~ 70 kms~^, or skewed 



towards higher velocities {^ 200 km s ) but still 



As discussed by Lipunov (1992) 



isolated NSs can 
be classified into four main types: Ejectors, Pro- 
pellers, Accretors and Georotators. In Ejectors 
the relativistic outflowing momentum flux is always 
larger than the ram pressure of the surrounding ma- 



rich in slow stars ( Paczyhski 1990| ). In the last few 
years however, proper motion studies and revised 
distance estimates of pulsar subpopulations revealed 
that young neutron stars have high mean velocities 
~ 400 km s 



by Lyne fc Lorimer (1994) 



^ according to the original suggestion 
200 km s~^ on the 



or ~ 



basis of the more recent analyses by Hansen & Phin- 



terial so they never accrete and are either active ra- ney (1997), Cordes fc Chernoff (1997) and Cordes & 



dio pulsars or dead pulsars, spun down by dipole Chernoff (1998). Although little is known about the 



losses. In Propellers the incoming matter can pen- 
etrate down to the Alfven radius but not further be- 
cause of the centrifugal barrier, and, although sta- 
tionary inflow can not occur, the piling up of the 
material at the Alfven radius may give rise to (sup- 
posedly short) episodes of accretion (Treves, Colpi 
& Lipunov 1993 ; Popov 1994| ). Steady accretion 
is also impossible in Georotators where (similarly to 
the E arth) the Alfven radius exceeds the accretion ra- 



low-velocity tail of the distribution (Hartman 1997; 
Hartman et al. 1997] ), the possibility that the low- 
velocity tail is underpopulated with respect to what 
was previously assumed should be seriously taken into 
account. It is our aim to revise the estimates on the 
number of old accreting neutron stars in the Galaxy, 
and in the solar vicinity in particular, in the light of 
these new data, in the attempt to reconcile theoretical 



dius, po that magnetic pressure dominates everywhere 



over the gravitational pull. It is the combination of 
the star period, magnetic field, velocity and ambient 
medium, that decides which type a given isolated NS 
belongs to and, since both P, B, V, and n change 
during the star evolution, a NS can go through dif- 
ferent stages in its lifetime. This argument shows 
that, at variance with what was assumed in earlier 



predictions with present ROSAT limits (Neiihauser & 
Triimper 1999| ). This is of potential importance as it 



invest igations (Treves fc Colpi 1991 ; Blaes fc Madau 
19931^ |Zane et al. 1995]) stationary accretion onto an 
isolated old NS (ONS hereafter) depends crucially on 
its rotational period and magnetic field. 

The dynamical evolution of NSs in the Galactic 
potential was studied by sever al authors ([Paczyhski 
1990|;]piaes fc Rajagopal 199 1|; [Blaes fc Madau 1991 ^ 



may indirectly reveal how the long term evolution of 
key parameters, such as the magnetic field, the pe- 
riod, the velocity distribution and the star formation 
rate relate to the properties of neutron stars at birth 
and to their interaction with the Galactic environ- 
ment. 

2. THE MODEL 

In this section we summarize the main hypoth- 
esis introduced to describe the evolution of single 
stars and outline shortly the technique used to ex- 



plore their statistical properties, referring to Popov 
[fc Prokhorov (1998)[ and Popov fc Prokhorov (1999 



Madau k Blaes 1994[ ; [Zane et al. 1995D who mod- 
elled the present velocity distribution, one of the key 
parameters governing accretion. Until now, how- 



for details on spatial evolution calculations. 

2.1. Dynamical evolution 

The dynamical evolution of each single star in the 
Galactic potential is followed solving its equations of 



Popov et al. 



motion. The potential, proposed by Miyamoto & Na- 
gai ( ^975) and adopted later by Paczynski (1990), is 
given by the superposition of a spherical halo and 
of two flattened bulge/disk components. In Galacto- 
centric cylindrical coordinates R, Z, r = ^E? + Z"^ 
they are expressed as 



^h{R.Z) 



GMc 

re 



1 / r' 
- In 1 + ^ 



r, 



and 



^i{R,Z) 



c / 
r \rc 



GMi 



i?2 + [a, + (Z2 + 62)l/2]2 



(1) 



(2) 



where the index i stands both for B (bulge) and D 
(disk) and H for halo. The values of the parameters 
appearing in equations (||) and (^) are summarized 
in table ^. 

In their motion through the Galaxy, NSs interact 
with the interstellar medium (ISM) when (and if) 
they enter the Propeller or the Accretor phase. In 
these two stages the accreting material affects signifi- 
cantly the star spin because braking torques arise (see 
§2.3). Since the period evolution depends on both the 
star velocity and the local density of the interstellar 
medium, any attempt to investigate the statistical 
properties of the NS population should incorporate a 
(detailed) model of the distribution of the ISM. Un- 
fortunately the distribution of molecular and atomic 
hydrogen in the Galaxy is highly inhomogeneous and 
only an averaged description is possible. Here we use 



the analytical distributions from Bochkarev (1992) 
and Zane et al. (1995) . Denoting with nui and n^/g 
the neutral and molecular hydrogen number density, 
the total proton density is given by 



n{R, Z) = nni + 2nH2- 



(3) 



The molecular hydrogen distribution is approximated 
as 

__^2 

UH^iR^Z) = n2{R) exp ^ _ ('jqj^^\2 



(4) 



where n2{R) is tabulated in Bochkarev (1992) . The 
map of atomic hydrogen is more complex and for 
R < 3.4 kpc we assume 



riHi = no (i?) exp 



-Z^ f R y 
2- (14000)2 UkpcJ 



(5) 



rotating ring and dense clouds. However, since we do 
not calculate the evolution of NSs born in the cen- 
tral 2 kpc, this is not going to be of any relevance 
Zane, Turolla fc Treves 1996a| for a discussion of 



see 



the diffuse X-ray emission from accreting NSs in the 
Galactic Center). The density of cold and warm HI 
at 3.4 kpc < R < 8.5 kpc can be fitted by 



0.107 exp 



nni = 0.345 exp 

-Z2 



2- (530pc)2 



2- (212pc)2 
-Z 



+ 



+ 0.064 exp 



403 pc 



.(6) 



while at radii larger than 8.5 kpc we use 



riHi = n3{R)exp 



-Z2 



2-(530pc)2 V8.5kpc 



R 



(7) 



where no{R) and n^^R) are again taken from 
Bochkarev (1992). The total hydrogen density dis- 



tribution used in our computations is shown in Fig- 
ure [l|; the density in the galactic plane varies by 
more than one order of magnitude, ranging from 
0.1 cm~^ to 4.3 cm~^. Within a region of ~ 140 pc 
around the Sun, the ISM is underdense, and we take 
0.07 cm"^ 



n 



(see 



Zane et al. 1996bD 



In our model we assume that the NS birthrate is 
constant in time and proportional in magnitude to 
the square of the local gas density (see for example 
Cox 19831 ; [Firmani fc Tutukov 199^ ). This imphes 



that NSs are preferentially produced in the molecu- 
lar ring region and migrate through the Galaxy dur- 
ing their evolution, due to the large kick velocities 
acquired at birth. Stars are assumed to be born 
in the Galactic plane {Z = 0) and in the range 
2 kpc < R < 16 kpc; this limitation has no ma- 
jor influence on the results, and on the number of 
sources found in the solar vicinity in particular. In 
the present investigation all effects of NSs born in 
binary systems (see, e.g., Iben fc Tutukov 1998| ) has 
been neglected. Mirroring the birthrate of massive 
stars, also the NS rate was likely to be higher in past 
( [White &: Ghosh 1998| ), so that our assumption of a 
constant NS birthrate is just a simplification. The 
effect of a time dependent NS formation rate will de- 
serve further study. 

Neutron stars at birth have a circular velocity de- 
termined by the Galactic potential. Superposed to 
this ordered motion a kick velocity V is imparted in 
a random direction. The exact form of the kick distri- 



bution is still uncertain (see, e.g., Lipunov, Postnov & 



The ;)receding expression becomes inaccurate inside 



Prokhorov 1996a; Iben & Tutukov 1998; Belczyhski 
fc Bulik 1999| ). Here we use an isotropic Gaussian dis 



1 kpc, where the gas distribution is dominated by a tribution with zero mean velocity vector relative to 
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the local circular speed, and one-dimensional veloc- 
ity dispersion ay, simply as a mean to model the true 
pulsar distribution at birth (see e.g. pordes 1998| ). 
The resulting distribution for the modulus V is char- 
acterized by a mean velocity (V) = (8/7r)"^'^(Ty which 
is varied in the interval 0-550 km s~^. The veloc- 
ity dispersion of the NS progenitors, ~ 20 km s~^, 
has been neglected and runs with (V) of this order 
mimic a NSs population with zero mean kick velocity 
at birth. 

For each star the dynamical evolution has been fol- 
lowed up to present time. The evolution of the spin 
period P (as described in §2.3) is also computed, 
keeping memory of the different phases the neutron 
star experiences (Ejector, Propeller, etc.). A number 
of representative evolutionary paths have been cal- 
culated, placing the stars at a given location in the 
galactic plane, with a randomly oriented kick drawn 
from the specified Gaussian distribution. Each evo- 
lutionary track is statistically independent, and gives 
a unique galactic orbit together with a unique his- 
tory. Each orbit is recorded every time the star 
crosses the boundary between two adjacent spatial 
cells in which the Galaxy is divided. The time of 
tion between two phases is also recorded and 



10 kms^^. Here and in the following M and R de- 
note the NS mass and radius, which we take equal to 
1.4 Mq and 10 km, respectively, for all stars. Assum- 
ing, for the sake of simplicity, a non-relativistic effi- 
ciency ~ GM/Rc^, the accretion luminosity is given 
by 



L 



GMM 
R 



(9) 



Although the spectrum emitted by accreting NSs may 
differ rather substantially from a blackbody at the 



star effective temperatu re (see Zampieri et al. 199£ 
for low-field spectra and|Zane, Turolla fc Treves 1999 



for magnetized models) here we assume pure thermal 
emission. Even in this simplified case, the effective 
temperature depends on the star dipole field because 
accretion occurs only at the two polar caps of radius 

Rr. 



''Cap 



R^{R/Ra), 



Reap = 9.5 X 10^ ^gj 



2/7.1/7,-3/7 



n 



-"10 



cm; 



(10) 



here Ra ^ 1.1 x 10 n 



-2/7„,6/7 4/7 



^10 



/igo cm is the Alfven 



trans: 



radius and //30 is the star magnetic dipole moment in 
units of lO^'' G cm^ (see e.g. Zane et al. 1996b; Ko 



store£ Approxhnatebr, a sti^ takes lO'^ 3T to cross a ^^"'f ^ ^, P°P°^ ^997|). The reduced emitting area 



cell, which means that the star position (and status) 
is recordered ~ 10 times in a total orbital time of 
"I^qIO yj. -jiq j^ijnig continuous star formation a single 
track is repeatedly used shifting the time ahead. In 
addition, tracks born in denser regions of the galactic 
plane, where the rate of star formation is higher, ac- 
quires, accordingly, a higher statistical weight. One 
orbit integration samples the behaviour of nearly 10^ 
stars born with the same initial conditions at differ- 
ent times, (we refer for details to the "Scenario Ma- 
chine" of [Lipunov, Postnov Sz Prokhorov 1996b ). A 
number of independent runs have been performed us- 
ing different number of objects, up to a maximum of 
~ 10^ so giving ~ 10^ stellar paths. The runs have 
been carried out using different seeds for the random 
orientation of kicks. Statistical fluctuations on the 
number of stars in each stage are typically of a few 
percent. 

2.2. Accretion physics 

The accretion rate was calculated according to the 
Bondi formula 



where nip is proton mass, the sound speed Vg is al- 
ways 10 kms~^ and vio = (V^ + V^)^^"^ in units of 



produces harder spectra, with an effective tempera- 
ture ~ 3-4 times larger than in the unmagnetized 
case 



r,;y=.5xl0V3iV3A%-39/i4K: 



(11) 



typical temperatures are around 100 eV for V ~ 100 
kms~^ for non-evolving magnetic field, and a little 
bit lower for evolving magnetic fields. 

2.3. Period evolution 

All neutron stars are assumed to be born with a 
period -P(O) = 0.02 s, and a magnetic moment ei- 
ther /i3o = 1 or /iao = 0.5. Different distributions of 
the initial periods, like the one recently proposed by 
[Spruit &: Phinney (1998) , were also tested and pro- 
duced very similar results. 

The ejector regime begins with the pulsar phase 
and proceeds also after the breakdown of the co- 
herence condition when the star becomes a dead, or 
silent, pulsar. In this phase the energy losses are due 
to magnetic dipole radiation and the period increases 
in time according to 



P 



P(0) + 3x IO-V30 t^^"^ s 



(12) 



where t is in yr. When the gravitational energy den- 
sity of the incoming interstellar gas exceeds the out- 
ward momentum flux at the accretion radius, Rac — 
2GM/v'^, matter starts to fall in. For this condition 
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to be met, the period must have reached a critical 
value 



P£;(^^P)~10^3^ n 



1/2 -1/4 1/2 



^10 



(13) 



which is attained by dipole braking (for a constant 
field) in a time 



tE ~ 10^ /i3o^ n 



-1/2 



^10 yr . 



(14) 



When P > Pe{E -^ P) matter can penetrate down 
to the Alfven radius, but the interaction with the ro- 
tating magnetosphere prevents accretion to go any 



further because of the centrifugal barrier. The NS is 
now m the propeller phase, rotational energy is lost 
to the ISM and the period keeps increasing at a rate 



dP 

'dt 



MR\P 



I 



KP'^ ss" 



(15) 



Here we take K = 2.4 x 10 



-14„8/7 3/7„-9/7 
""" " "lO 



a 



1 



(Shakura 1975) and / = 10'^^ gcm^ is the star mO' 



Ho n' 
gem 

ment of inertia. 

It should be stressed, however, that our expres- 
sion for the propeller spin-down is just an approx- 
imation. The propeller physics is very complicated 
and its thorough understanding requires a full 2-D 
or even 3-D MHD investigation of accretion onto a 
rotating dipole (see e.g. Toropin et al. 1999t see also 
other spin-down formulae for that stage in Lipunov 
& Popov 1995). The propeller spin-down, as mod- 
eled by equation (|l^), is very efficient and acts on a 
typical timescale 



tp ~ 1.3 X 10°;U3o 



8/7„-3/7„9/7 



n 



10 yr 



(16) 



Numerical simulations ( [Toropin et al. 199£ ; Toropin, 
private communication) seem indeed to confirm that 
spin-down in the propeller phase is very fast due to 
the large mass expulsion rate. Note that the present 
expression for the propeller torque is somewhat differ- 
ent from that adopted by Treves et al. (1993; 1998). 
As the star moves through the inhomogeneous ISM a 
transition from the propeller back to the ejector phase 
may occur if the period attains the critical value 



PE{P^E)^3fsl', 



4/5 6/7 -2/7 



(17) 



Note that the transitions P ^ E and E ^ P are not 
symmetric because the two periods (|^) and (17) de- 
rive from the same physical condition but evaluated 
at different radii, as first discussed by Shvartsman in 
the early '70s. 

Accretion onto the star surface occurs when the 
corotation radius Rco = {GM P'^/Att'^)^'^ becomes 
larger than the Alfven radius (and Ra < Rac-, see 



below). This implies that braking torques have in- 
creased the period up to 



PA(i^^vl) -420^3^ n 



6/7„-3/7„9/7 



^10 



(18) 



As soon as the NS enters the accretor phase, 
torques produced by stochastic angular momentum 
exchanges in the ISM slow down the star rotation to 
the equilibrium period 



P. 



eq 



2.6 X 10^ t; 



,3 ,,-2/3, 2/3 ^-2/3 ,,13/3 



(t)10 



^30 



n 



^10 



(19) 



where f (^^ the t urbulent velocity of the ISM (Lipunov 
Popov 1995] [Konenkov &: Popov 19971 ). 



Actually the condition that P > Pa is not sufficient 
to guarantee that matter is captured at the accretion 
radius. At the very low accretion rates expected for 
fast, isolated NSs, it could be that the Alfven radius 
is larger than the accretion radius. The condition 
Ra < Rac translates into a limit for the star velocity 



7;<410n^/^V3o^^^ kms" 



(20) 



Above this velocity the NS behaves as a georotator 
and it experiences a torque which is computed in the 
same way as in the propeller stage. The star can en- 
ter the georotator phase only from the propeller or 
accretor stage. 

Figure |2| illustrates the various possible stages for a 
NS lO^'' yr old as a function of the star velocity and 
magnetic field for a constant ambient density of 1 
cm-3. The propeller region in Figure |2| is very small, 
as a consequence of the extremely efficient spin-down 
implied by equation (ll5|) with our choice of the pa- 
rameters Q and K. 

3. RESULTS AND DISCUSSION 

In this section, we present the results of our numer- 
ical simulations both in the cases of a constant and 
of a decaying magnetic field. 

3.1. The NS census for a non-decaying field 

Here we consider two values for the magnetic dipole 
moment (//30 = 0.5 and /U30 = 1) as representative of 
the constant NS magnetic field. The present frac- 
tion of NSs in the Ejector and Accretor stages as a 
function of the mean kick velocity is shown in Fig- 
ure y. Statistical errors are typically ~ 2%. Pro- 
pellers and Georotators are much less abundant and 
never exceed ~ 1% of the total number. In addition, 
their fractions oscillate intrinsically because of the 
changes in the velocity and in the ISM density along 
the star trajectory. These effects add further noise 
to the statistical fluctuations and, given the small 
number of sources in these two states, prevent any 
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definite conclusion about the dependence of the frac- 
tions from the mean kick velocity. We note that the 
fraction of Accretors increases with the field strength, 



in ag'eement with the findings of Livio, Xu & Frank 



unusually high initial magnetic fields or/and long pe- 
riods ( potthelf 1998 ), or to the fall-back in the af- 
termath of the supernova explosion (|Colpi, Shapirc] 



(1998|)| and |Colpi et al. (1998)| . Highly magnetized 



Wassermann 1996; Gepper, Page & Zannias 1999) 



NSs suffer, in fact, a more severe spin-down during 
the Ejector and Propeller phases and reach the accre- 
tion phase earlier in their history. For large enough 
mean kick velocities (> 400 km s^^ for fi^Q = 1 and 
> 300 km s~^ for fi^Q = 0.5) the fraction of accretors 
becomes very small and again results become statis- 
tically uncertain. We note also that the numerical 
simulation produces a picture of the present status of 
the Galactic NS population which is consistent with 
that emerging from Figure |2|. In fact, if we assume 
that the velocity gives a measure of the average kick, 
we see, with reference to /iso = 1, that Accretors are 
more abundant than Ejectors up to '^ 100 kms~ , 
Propellers are extremely rare and Georotators nearly 
absent. This is precisely what Figure |3| shows. 

To obtain an estimate of the number of accreting 
sources in the Solar vicinity (taken to be a sphere of 
radius 140 pc centered on the Sun), we used all the 
stars contained in a torus of the same radius located 
at 8 kpc from the center of the Galaxy. The result 
was then rescaled to the volume of interest and is 
shown in Figure 0. Here, and in the following the 



Radio pulsars are observed only in a fraction of SNRs 
(see e.g. Kaspi 1996 ; Frail 199q ), and even some of 
these coincidences are doubtful, so there is a serious 
possibility that the total number of NSs derived from 
radio pulsar statistics is only a lower limit. 

In order to compare the expected number of accret- 
ing ONSs with the ROSAT Ah Sky Survey (RASS) 
results, we evaluated the number of those ONSs, 
within 140 pc from the Sun, producing an unabsorbed 
flux of 10"^^ erg cm~^ s~^ at energies ~ 100 eV (the 
contribution from sources at larger distances is found 
to be negligible). This flux limit should correspond 
to a RASS count rate of 0.01 cts s~^ for a black- 
body spectrum, once the interstellar absorption and 
the ROSAT response function are accounted for. The 
results are illustrated in Figure ^ (above 200 km s~^ 
statistical errors are dominant). As expected, the 
dependence of the number of visible sources on the 
kick velocity is rather strong. The main result is that 
for mean velocities below 200 km s~^ the number 
of ONSs with a flux above the RASS detection limit 
would exceed 10. Most recent analyses on the number 



of isolated NSs in the RASS ([Neiihauser fc Triimper 



total number of Galactic NSs was assumed to be 10^. 1999| ) indicate that the upper limit is below 10. This 



As expected, the local density is a sensitive function 
of the kick velocity and varies between n ~ 8.5 x 10~^ 
pc^'^ and n ~ 6 X 10~^ pc~^ for (V) = and 190 
km s~^ respectively. These figures should be com- 
pared with n ~ 1.4 X 10"'^ (A^/10^) pc~^ as derived 
by Paczyiiski (1990) for a zero mean velocity and with 
n ~ 3 - 7 X 10'"^ (N/W^) found by Blaes & Madau 
(1993) and Zane et al. (1995) for an initial mean ve- 
locity of ~ 60 km s~^. The values of the local density 
of isolated NSs are thus consistent with these previ- 
ous estimates. 

The density of isolated NSs projected onto the 
Galactic plane is ~ 2.4 x 10^ (A^/10^) kpc~^ for aver- 
age initial velocities ~ 200 km s~^ and a scaleheight 
of 200 pc. This figure is close to the one deduced 
by Neiihauser &: Triimper (1999) from radio pulsars 
observations (Lyne et al. 1998) with the assumption 
of a pulsar lifetime ~ 10^ yr. We would like to note, 
however, that a projected density ~ 10^ kpc^^ is now 
obtained for N = 1 0^ instead of 10^ (as in Madau & 



Blaes 1994 ; see also NFeiihauser &: Triimper 1999 ) be- 
cause of the larger mean velocity at birth. A total 
number ~ 10^ appears to be consistent with the nu- 
cleosynthesis and chemical evolution of the Galaxy, 
while 10^ is derived from radio pulsars observations. 
For the time being, it is still uncertain if all NSs ex- 
perience an active radio pulsar phase, due to low or 



implies a lower limit for the mean kick velocity at 
birth of ~ 200 km s~^ for a total number of stars 
~109. 

This result is in agreement with the estimates de- 
rived from pulsar statistics. An important aspect is 
that our results exclude the possible presence of a 
low-velocity tailat birth in excees to that contained 
in a Gaussian with {V) ~ 200 km s~^, i.e., the 0.05% 
(= (27r)~(^/^)[F/(T]'^) of the whole population, when 
considering stars having U < 40 km s^^. 

3.2. The NS census for a decaying field 

The time evolution of the magnetic field in isolated 
NSs is still a very controversial issue and no firm con- 
clusion has been established as yet ([Bhattacharya et 
al. 1992| ). A strong point is that radio pulsar obser- 
vations (see e.g. Lyne et al. 1998 ) seem to rule out 
fast decay with typical times less than ~ 10 Myr, but 
this does not exclude the possibility that B decays 
over much longer timescales {td ~ 10^ — lO^'' yr). For 
this reason we have investigated to what extent the 
decay of the S-field influences the results presented 
in the previous section. We refer here only to a very 
simplified picture in which B{t) = B{0)exp {—t/td) 
and no attempt is made to justify this law on a phys- 
ical ground. (See Popov &: Prokhorov 1999] for a more 
complete study in the parameter space.) In this re- 
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spect we just mention that detailed models predict 
both exponential (quite similar to the one assumed 
here) and non-exponential decay (e.g. a power-law, 
Urpin & Konenkov 1997). Calculations have been 



performed for t^ = 1.1 x 10 yr, td = 2.2 x 10 yr and 
^30 (0) = 1. Since no bottom field was specified, the 
magnetic moment becomes very low for stars born 
~ 10^° yr ago (~ 10^^ G cm^ and ~ 10^^ G cm^ for 
the two values of td respectively) . 

Results are summarized in Figure p. As it is ex- 
pected (see Colpi et al. 1998), the number of Pro- 
pellers is significantly increased with respect to the 
non-decaying case, while Ejectors are now less abun- 
dant. Georotators are still very rare, ^ 1%, and are 
not shown in Figure IBI. The fraction of Accretors is 
approximately the same for the two values of td, and, 
at least for low mean velocities, is comparable to that 
of the non-decaying field while, at larger speeds, it 
seems to be somehow higher. This shows that the 
fraction of Accretors depends to some extent on how 
the magnetic field decays. We know that a core field 
decaying over a time ~ 10*^ yr before freezing would 
produce an underabundance of Accretors relative to 
the case of a constant field (Colpi et al. 1998) be- 
cause NSs persist in the Ejector or Propeller phase 
never approaching the Accretor phase. By contrast, a 
fast and progressive decay of B would lead to an over- 
abundance of Accretors because this situation is simi- 
lar to "turning off" the magnetic field, i.e., quenching 
any magnetospheric effect on the infalling matter. 

Summarizing, we can conclude that, although both 
the initial distribution and the subsequent evolution 
of the magnetic field strongly influences the NS cen- 
sus altering the fraction of Ejectors relative to Pro- 
pellers, the lower bound on the average kick derived 
from ROSAT surveys is not very sensitive to B, at 
least for not too extreme values of td and ^(0), within 
this model. 



for a constant or a decaying 5-field, unless the decay 
timescale is shorter than ~ 10^ yr. Since observable 
accretion-powered ONSs are slow objects, our results 
exclude also the possibility that the present velocity 
distribution of NSs is richer in low- velocity objects 
with respect to a Maxwellian. The paucity of ac- 
creting ONSs seem therefore to lend further support 
in favor of neutron stars as being very fast objects 
spending most of their live in the Ejector phase. 
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4. CONCLUSIONS 

In this paper we have investigated how the present 
distribution of neutron stars in the different stages 
(Ejector, Propeller, Accretor and Georotator) de- 
pends on the star mean velocity at birth. On the 
basis of a total of ~ 10^ NSs, the fraction of Ac- 
cretors was used to estimate the number of sources 
within 140 pc from the Sun which should have been 
detected by ROSAT. Most recent analysis of ROSAT 
data indicate that no more than r^ 10 non-optically 
identified sources can be accreting ONSs. This im- 
plies that the average velocity of the NS population 
at birth has to exceed ~ 200 kms~^, a Figure which 
is consistent with those derived from radio pulsars 
statistics. We have found that this lower limit on the 
mean kick velocity is substantially the same either 
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Table 1 
Parameters for the various contributions to the Galactic potential 

Halo Bulge Disk 

re = 6.0 kpc qb = 0.0 kpc ao = 3.7 kpc 

bB = 277 pc bn = 200 pc 
Mc = 5.0 X 10^° Mq Mb = 1.1 x 10^° M© Md = 8.1 x 10^° Mq 
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Fig. 1. — The hydrogen number density distribution in the R-Z plane. 
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Fig. 2. — The different stages of old NSs at present as a function of the star velocity and magnetic dipole moment. 
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Fig. 3. — Fractions of NSs in the different stages vs. the mean kick velocity for /iso — 0.5 (triangles) and /130 — 1 (diamonds); 
typical statistical uncertainty is ~ 2%. 
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Fig. 4. — Total number of NSs in the Solar vicinity {R < 140 pc) for fi,io = 1. 
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Fig. 5. — Number of accreting NSs in the Solar vicinity above the ROSAT All Sky Survey detection limit for a constant (/lao = 1, 
diamonds) and decaying field {td — 2.2 x 10^ yrs, triangles). 
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Fig. 6. — Fractions of NSs in the diflerent stages vs. the average kick velocity for a decaying field with an e-folding time 
td — 2.2 X 10^ yrs (triangles) and td — 1.1 x lO'"* yrs (diamonds). 



